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Radiative decays of X(3872) are studied in single-channel approximation (SCA) and in the 
coupled-channel (CC) approach, where the decay channels DD* are described with the string 
breaking mechanism. In SCA the transition rate IY = F(2 3 Pi —> 0y) = 71.8 keV and large 
fi = r(2 3 Pi -4 J/ip'y) = 85.4 keV are obtained, giving for their ratio the value R^~, = & = 0.84. 

In the CC approach three factors are shown to be equally important. First, the admixture of 
the l 3 Pi component in the normalized wave function of X(3872) due to the CC effects. Its 
weight cx(Pr) = 0.200 ± 0.015 is calculated. Secondly, the use of the multipole function g(r) 
instead of r in the overlap integrals, determining the partial widths. Thirdly, the choice of 
the gluon-exchange interaction for A (3872), as well as for other states above threshold. If for 
A'(3872) the gluon-exchange potential is taken the same as for low-lying charmonium states, then 
in the CC approach Ti = F(A(3872) —> J/ip'y) ~ 3 keV is very small, giving the large ratio 
R^, 1 = 10. Arguments are presented why the gluon-exchange interaction may 

be suppressed for A(3872) and in this case Fi = 42.7 keV, lA = 70.5 keV, and P^, 7 = 1.65 are pre¬ 
dicted for the minimal value cx(min) = 0.185, while for the maximal value cx = 0.215 we obtained 
Ti = 30.8 keV, T 2 = 73.2 keV, and R/, 7 = 2.38, which agrees with the LHCb data. 

PACS numbers: 11.10.St, 12.39.Pn, 12.40.Yx 


I. INTRODUCTION 

In 2003 the Belle collaboration discovered the A(3872) as a narrow peak in the J/ipirn invariant mass distribution in 
the decays B -4 J/ipTntK [l]]. Now its characteristics, like the mass, the strict restriction on the width, T < 1.2 MeV, 
and the charge parity C = +, are well established In recent CDF and LHCb experiments the quantum 

numbers of A(3872) were determined to be J PC = 1 ++ 10. Still, discussions about the nature of X(3872) 
continue and to understand its exotic properties a special role is played by the radiative decays, X(3872) -4 J/ip'y 
and X(3872) -4 0(3686)7, which are sensitive to the behavior of the A'(3872) wave function (w.f.) at medium and 
large distances. 

The first evidence for the decay X(3872) -4 J/ip'y was obtained by the Belle collaboration Q and confirmed by the 
BaBar collaboration 0. Later BaBar has also observed the radiative decay X(3872) -4 0(3686)7 and determined the 
branching ratio fraction R^, 7 = B ^(x( 3872 )^f/^ 7 ) 7 ^ = 3-4 ± 1.4 [l0j. However, Belle has not found evidence for the 
radiative decay X(3872) -4 ■0(3686)7 and put an upper limit for the ratio & 

R in < 2.1. (1) 

Recently the LHCb group [l2j has observed the decay X(3872) -4 0(3686)7 w ith a good statistics and determined 
its value to be 


= 2.46 ± 0.64 ± 0.29, (2) 

while in Ref. fl3} the weighted average over three groups of measurements was determined to be R^ )1 = 2.31 ± 0.57. 

Unfortunately, theoretical predictions for the partial widths Ti and T 2 of the radiative decays X(3872) -4 J/ip 'y 
and X(3872) -4 0(3686)7, respectively, vary widely in different models flil-lxil (see also the recent reviews [23, 241]). 
If X(3872) is considered as a pure 2 3 Pi charmonium state, then a large value R^ ~ 5 is obtained as shown in 
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Refs. [l6l4l9j . exceeding the LHCb result. On the contrary, in a molecular picture the transition X(3872) — > if>( 3686)7 
is suppressed and the ratio should be much smaller urn!- Also in many-channel models, both in the 3 Pq 
model ill, ij and the Cornell model [f9|, large values i ?^ 7 > 5.0 were predicted. 

It is surprising that the predicted widths lb and T 2 strongly differ even within the single-channel approximation 
(SCA), i.e., when X (3872) is supposed to be the 2 3 Pi charmonium state [TgI - ITqI . It occurs because different parameters 
in the potentials and different kinematics are used. Therefore it remains unclear to what degree the theoretical 
predictions refer to specific features of the radiative decays or to the choice of fitting parameters and kinematics. 
Meanwhile, during the last decade within the field correlator method @ - (27j and in lattice QCD @, @| the static 
potential was shown to be universal in the region r < 1.5 fm, where there is no creation of light qq pairs. Here we shall 
use this information and perform a parameter-free analysis of the charmonium states with the use of the relativistic 
string Hamiltonian (RSH) [25j |. 

It was already underlined by Li and Chao and in Ref. |24j that the use of relativistic kinematics is very important 
for the radiative decays and in nonrelativistic approximation the overlap integrals, determining the partial widths, 
may significantly differ from those in relativistic calculations. 

In our previous paper 3 the rates and the ratio i ?^ 7 were calculated considering X(3872) (with J PC = 1 ++ ) 
in the coupled-channels (C C) app roach, where the coupling to the DD* channels is determined by the relativistic 
string-breaking mechanism |20r-22j. This mechanism was successfully applied to X(3872), explaining it as the 2 3 Pi 
charmonium state shifted down due to CC effects and appearing as a sharp peak just at the D°D*° threshold |20j | . 
Owing to such CC effects the w.f. of X (3872) acquires an admixture cx(Pr) from the 1 3 Pi component, equal to 
~ 0.20. Although this mixing parameter is not very large, it nevertheless strongly affects the value of the partial 
width Ti, in which the overlap integral is small and very sensitive to different corrections. 

However, in Ref. @] in the overlap integral Ki = (ip(2P)\g(r)\ip(lS)) the multipole moment g(r) was replaced by 
r, which is a common practice. In this approximation A'i = (y>(2P)|r’|y>(lS)) = 0.21 has a small value, close to those 
given in Refs. [H|,3- It appears that the magnitude of K\. defined using the multipole function g(r), is a factor of 
two larger than K 1 , increasing the partial width Ti by about four times. Thus the replacement of g(r) by r cannot 
be applied to this radiative decay, proceeding with a large photon energy. 

We also show that the integral K\ strongly depends on the gluon-exchange (GE) part of the static potential and 
consider two possibilities: First, when in A(3872) the GE potential is taken the same as for the low-lying states (below 
threshold). In this case the results are different in the SC A and in the CC approach. The second possibility implies 
that in higher charmonium states (above threshold), including A(3872), the GE interaction is suppressed. For such 
a cc dynamics the partial width Fi increases, while the partial width T 2 = T(A(3872) —> -^( 3686 ) 7 ) weakly depends 
on the GE interaction. 

We use here the CC approach, where the admixture from the continuum to the w.f. of X (3872) is defined by the 
mixing parameter cx(Er = 3.872 GeV), which is calculated here with a good accuracy: cx = 0.200 ±0.015. Although 
its value is not large, it nevertheless gives an important contribution to Ti. 

Also in our analysis we lay the stress on the choice of the c-quark mass, the QCD constant A(n/ = 3), and the 
string tension, which cannot be arbitrary parameters lest the physical picture would be distorted. 


II. THE cc INTERACTION BELOW AND ABOVE THE THRESHOLD 


Before studing the X(3872) in the many-channel approach one needs to know its w.f. in the SCA. The rate of the 
electric dipole (El) transition between an initial state i: n 3 P\ and a final state /: m 3 S\ is given by the expression 

@13, 


T(i^f + 1 ) = i a e 2 Q k 3 (2J f + 1) Sf f \£ if \ 2 , (3) 

which contains the overlap integral £)/ (called also the m.e.), calculated here via the relativistic w.f. The photon 

M?—M 2 

energy fc 7 is defined as fc 7 = — ' 2M f , where Mj(Mj) is the mass of the initial (final) state. In Eq. ([3]) the statistical 
factor Sf ; = S®- is given by 


Sf f = max (l, V) 


J 1 J' \ 

Vs l J 5 


(4) 


which for the El transitions between the n 3 Pj and m 3 Si (m 3 Di) states, with the same spin S = 1, is equal to 

I (A) 

9 vi8v 
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The overlap integral /, given by 

£ if = J drr 2 R ni pr)g{r)R nf j t (r) = {n i ,l i \g{r)\n f ,l f ), 
is determined by the relativistic radial w.f. R ni ,ii{f) and R nfi i f (r). The function g(r) [ 


is: 


g{r) = T~[yjo{y) ~ ji(y)\ = “ {sin y (1 - \ ) + - cosy 

2 y { V r / y 


( 5 ) 


( 6 ) 


where j n (y) (n = 1, 2) are the spherical Bessel functions and the variable y = kjr/2. It can be easily shown that the 
replacement of g(i') by r ( g(r ) < r), used in many analyses, is a good approximation for y < 1.0, while in the range, 
1.1 < y < 1.3, the difference g(r) — r can already reach < 30%. The zero of the function g(r) occurs at z/ 0 = 2.74, 
corresponding to very large r, r > 10 GeV -1 , even for the photon energy, fc 7 = 0.697 GeV for the V(3872) — > J/ip'y 
radiative decay. In r-space the difference g(r) — r remains small for photon energies fc 7 < 0.40 GeV. 

The overlap integrals, involving the IP and 2 P states, are denoted below as 


Im — IlP,mS — 

J drr 2 Rip{r)g(r)Rm.s(r), 

(m = 1,2), 


K-m = I^-2P,mS = 

1 drr 2 R 2P (r)g(r)Rms(r), 

(to = 1,2). 

(7) 


In Eq. 0 the relativistic radial w.f. are calculated with the use of the RSH Ho, derived in Ref. [H] and applied many 
times to different effects in quarkonia nnm. For heavy quarkonia Hq has a very simple form: 


H a = 2^pi 


V B (r), 


( 8 ) 


where its kinetic term is similar to that in the relativized quark model (RQM) f32|. However, in Hq, by derivation, 
the fundamental value of the c-quark mass is equal to its pole mass, m c ~ 1.42 GeV. (It takes into account corrections 
perturbative in a s (m c ) and corresponds to the conventional current mass to c (to c ) ~ 1.23 GeV H3). In the RQM the 
c-quark mass is considered as a fitting parameter, e.g., m c = 1.628 GeV is used in Ref. [32] and m c = 1.562 GeV in 

Ref. [13- 

In the RSH the universal static potential V B (r) contains the confining and GE terms, 


V B (r) = ar - 


4q:b (r) 
3 r 


(9) 


where the string tension er = 0.18 GeV 2 is fixed by the slope of the Regge trajectories for light mesons and therefore 
cannot be considered as a fitting parameter. 

An important step is to make the correct choice of the vector coupling a B (r) j34|. To this end, the parameters 
defining the coupling, are here taken in correspondence to the existing data for the strong coupling a s (r ) in the 
MS scheme in perturbative QCD j35|. The asymptotic-freedom behavior of a B is determined by the “vector” QCD 
constant Ab, which, however, is not a new parameter but expressed through Rjjg- In particular, for n/ = 3 the 
relation Ab(u/ = 3) = 1.4753 A jjgijif = 3) is valid. In pQCD 


A Mg(n f = 3) = 339 ± 10 MeV, 


( 10 ) 


was extracted with the use of the matching procedure at the quark-mass thresholds [35|. This value Ajjg(nf = 3) = 
(339 ± 10) MeV gives a rather large value of Ab (rif = 3) = (500 ± 15) MeV, while in Ref. 341 from the analysis of the 
bottomonium spectrum, a smaller Ab (n/ = 3) = (465 ± 20) MeV was shown to be preferable; its value corresponds 
to A jig{nf = 3) = (315 ± 14) MeV, in agreement with the lower limit in Eq. (flOl) . 

Also in Ref. [HJ] the infrared regulator M B was introduced, since for n/ = 3 the characteristic momenta q 2 are 
rather small and nonperturbative effects become important. As shown in Ref. [36|, the regulator M B is not a new 
parameter, but expressed via the string tension a, according to the relation Mg = 27rer, derived with an accuracy 
~ 10%. Then for a = 0.18 GeV 2 one finds that Mb = 1.06 ± 0.10 GeV; here the value Mb = 1.10 GeV is used. 

Thus the potential V B (r) contains only fundamental quantities: the current (pole) c-quark mass, the string tension 
from the Regge trajectories, and A = 3). In our calculations the following set of parameters is used: 

m c = 1.425 GeV, a = 0.18 GeV 2 , M B = 1.10 GeV, A B (n/= 3) = 465 MeV. 

For these values of As(n/ = 3) and Mb we find the frozen (asymptotic) value of the coupling a cr i t = 0.6086. 


( 11 ) 
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A remarkable property of the RSH is that in heavy quarkonium the centroid masses M cog (nl) of the nl multiplet 
exactly coincide with the eigenvalues of the spinless Salpeter equation (SSE): 

H 0 ip n i = M 0 (nl)<p n i. (12) 

In the charmonium spectrum not many levels are below the open charm threshold, where the universal static potential 
Vs(r) may be tested. However, above the DD threshold the cc interaction may differ from that for low-lying states, 
due to open channels or light-quark pair creation. In the string-breaking picture it is manifested as a flattening of 
the static potential at relatively large distances and this phenomenon is seen in lattice QCD [37j and also observed 
in the radial excitations of light mesons (38}. For that reason for higher states (above threshold) we consider two 
possibilities: the first one, when the GE potential is taken the same as in Eq. (0 and the second case, when the 
GE interaction is supposed to be suppressed for higher states, so that the cc dynamics is totally determined by the 
confining potential. There are several arguments in favor of the latter assumption: 

1. Suppression of the GE interaction follows from the analysis of the orbital Regge trajectories in charmonium 
[39|, |40[, which are not compatible with a significant GE contribution to the masses of higher states. 

2. The fine-structure effects, caused by the GE interaction, are not seen for the 2 P charmonium multiplet. For 
example, the measured mass difference 8M(2P) = M(y C 2(2P)) — M(y c o(2P)) < 10 MeV, is very small compared 
to the same mass difference for the IP multiplet: SM(1P) = M (x c 2 (lP)) — M (Xco(lP)) = 142 MeV [33], being 
smaller by an order of magnitude. It is difficult to explain such a suppression of the fine-structure effects, if a 
universal GE potential is used both for high- and low-lying states. (The suppression of the fine-structure effects 
does not contradict possible hyperfine splittings of the higher S— wave charmonium states due to a short-range 
hyperfine potential [271]). 

3. Our calculations of the charmonium spectrum show that the masses of higher states are practically the same 
for the linear and linear+GE potentials, if the correct choice of the c-quark mass is made (see below Table IJ). 

The confining potential used 

Vc(r) = or, (13) 

has the same string tension a = 0.18 GeV 2 . However, for the linear potential the c-quark mass may be different, since 
for a totally suppressed GE interaction there are no perturbative corrections to the c-quark mass and in this case m c 
is equal to the current mass (here m c (fh c ) = (1.27 — 1.29) GeV is used). 

In Table]!] we give the centroid masses M cog (nl) for both potentials, Eqs. (]7|) and (]9ji. 


TABLE I: The centroid masses M cog (nl) (in MeV) of higher multiplets for the potentials Vb (r) and Vc(r) 


State 

Potential Vb(t) 

Potential Vc(r) 

experiment 


m c = 1.425 GeV 

m c = 1.290 GeV 


3S 

4092 

4112 

4039 ± 1 

4S 

4447 

4448 

4424 ± 4 

2P 

3949 

3949 

3927 ± 3 a) 

3P 

4319 

4298 

absent 

ID 

3802 

3788 

3773 ± 3 

2D 

4187 

4155 

4153 ± 3 


For M cog (2P) the experimental mass M(\c2(2 3 P 2 ) is given. 


From Table]!] one can see that both spectra coincide within < 20 MeV accuracy. Moreover, the masses M(1D), 
M(2D), M(2P), and M(4S), calculated in SCA, are close to the experimental masses of 0(3770), 0(4160), y C 2 ( 2 P), 
and 0(4420), respectively, and one can expect that their hadronic shifts are not large, < 20 MeV. It is not so for the 
3 S state, where M(3S) is about (60 ± 10) MeV larger than M(0(4O4O)), indicating the importance of the CC effects 
for 0(4040). 

Thus our analysis of the charmonium spectrum does not allow to draw a definite conclusion what type of the cc 
interaction is preferable above threshold. Therefore other phenomena in charmonium have to be studied. One of them 
is just the radiative decay .¥(3872) —> J/0 7 , which, as was already underlined by E. Swanson m, “is particularly 
sensitive to model details”. 
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III. THE RADIATIVE DECAYS OF A(3872) IN SINGLE-CHANNEL APPROXIMATION 

Here we calculate the partial widths of the radiative transitions, A'(2 3 Pl) —> J/ip'y and AT(2 3 Pi) —> ip{2S)'y in 
SCA, i.e., assuming that A(3872) is the 2 3 Pl state. Then the widths (in units GeV) can be written as 


T? = 1.4418 10” 3 |A'i| 2 , (14) 

T° = 1.4418 10" 3 kl, |P' 2 | 2 , 

where the number 1.4418 10" 3 is the product of known factors in Eq. (O, £q 7 = 0.6974 GeV, k 2l = 0.1815 GeV, and 
for AT(3872) as the 2 3 Pl charmonium state, the integrals K\, K 2 (in the units GeV -1 ) are defined in Eq. ©■ 

From Eq. 0 and for the potential Vb(?') we obtained 

ATi =K 2Pt is = -0.418 GeV -1 , (15) 

which in SCA gives a large partial width: 


Ti = 85.4 keV, (16) 

This integral Ad, defined with g(r), has a magnitude that is two times larger than Ad = ( R 2 p\r\Ris), where g(r) is 
approximated by r: 

Ad = -0.210 GeV" 1 , id = 21.6 keV. (17) 

Notice, that even smaller magnitudes of K\ were obtained in other models [l7j], [21] (see Table ITT1). (The choice of 
the negative sign of Ad will be explained in Sect. ©. Thus, in the overlap integral determining the radiative decay 
X(2 3 Pi) — > J/ip 7 with rather large photon energy, fc 7 = 0.697 GeV, the replacement of g(r) by r has a poor accuracy. 
On the contrary, the differences between the overlap integral K 2 and the width T 2 , 

Ad = (2P\g(r)\2S) = 2.886 GeV" 1 , T 2 = 71.8 keV, (18) 

and the approximate quantities K 2 and T 2 

Ad = (2P|r|25) = 3.017 GeV" 1 , f 2 = 78.5 keV (19) 

are rather small, < 8 %. Then, because of the small value of id, the ratio of the partial widths T 2 to id is large, 


P^ 7 = 3.6. (20) 

On the contrary, if the integrals Ad and K 2 are calculated with g(r ): Ti = 85.4 keV, T 2 = 71.8 keV, then 

R 0 7 = 0.84, (21) 

is four times smaller. This ratio is not compatible with the BaBar [Oj and the LHCb results, Eqs. 0, but does not 
contradict the Belle restriction, P^, 7 <2.1. 

In Table mi we give the partial widths Ti and T 2 , predicted in several relativistic models. One can see that in the 
SCA, even with relativistic kinematics, the predicted values of Ti vary in a wide range, (11 — 85) keV. 

Surprisingly, the m.e. Ad in Eq. (1151) and Ti in Eq. (1161) calculated here with a relativistic Hamiltonian, practically 
coincide with the approximate m.e. ( 2P\r\\S) = 0.416 GeV" 1 and Ti = 84.6 keV, calculated in a nonrelativistic 
model with the logarithmic potential in Ref. [42]; the predicted value found there i?^ 7 = 1.1 is also close to our 
number in Eq. m- 


As a test of the w.f., defined by the RSH and used here, we also calculated the partial widths of the El transitions: 
\cj(lP) —> J/ip 7 (J = 0,1,2), where all three states lie far below the DD threshold and are described by the 
linear+GE potential. For % c i(lP) —> J/ip'y and Xc 2 (l-P) -A J/ip'y the calculated widths appeared to be in precise 
agreement with the experimental data. For these radiative decays the calculated integral hp,is = J(1 3 Pr./, Iff) = 
1.849 GeV " 1 is taken the same for J = 0,1, 2. For -0(3686) its w.f. contains an admixture from the ID state (43[ and 
therefore this radiative decay cannot be used as a test. 

For all three radiative transitions the photon energies are not large, fc 7 < 0.40 GeV, and therefore the differences 
between the m.e. (g(r)) and (r) are rather small, < 7%. Nevertheless, with the use of the function g(r) the partial 
widths of those radiative decays decrease by ~ 10%, improving the agreement with experiment. 

From Table mu one can see that the calculated values T(x c 0 (lP) —> J/ip'y) = 137 keV, T(x cl (lP) — > J/ip'y) = 
290 keV and T(x c2 (1P) — > J/ip'y) = 390 keV are in precise agreement with the experimental widths: r exp (x c o 
J/ip'y) = 133 ± 10 keV, T exp (xci -A J/ip'y) = 296 ± 13 keV, and T exp (x C 2 -5- J/ip'y) = 386 ± 17 keV (H- 
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TABLE II: The partial widths of the A'(3872) radiative decays (in keV) in relativistic models 


Transition 

fc 7 

[41] 

[16] 

[17] 

[18] 

our paper 

our paper 

Method 

in GeV 

SCA 

CCM 

SCA 

SCA 

SCA 

CCM 

2 3 Pi l 3 ^ + 7 

0.697 

33 

45 

11 

11 

85.4 

37 ±6 

2 3 Pi ->■ 2 3 Si + 7 

0.181 

146 

60 

70 

63.9 

71.8 

72 ± 1 

2 3 Pi -> 1 3 Di +7 

0. 098 

7.9 


4.0 

3.7 

5.8 



“1 The partial widths are calculated either in single-channel approximation (SCA) or within a coupled-channels method 

(CCM). 


TABLE III: The overlap integrals and the partial widths of the radiative decays y c j(lP) J/il’I- 


Transition 

fc 7 (GeV) 

(r) (GeV" 1 ) 

( 9 {r)) a > (GeV" 1 ) 

r (th) 

T(exp) 

1 3 P c 2 -t J/V’ + 7 

0.429 

1.976 

1.849 

390 

386 ± 17 

1 3 Pd -t J/ip + 7 

0.389 

1.976 

1.849 

290 

296 ± 13 

1 3 P c o —> J/ip + 7 

0.303 

1.976 

1.849 

137 

133 ± 10 


a ^ Calculations refer to the variant with m c = 1.425 GeV, Ab (n/ = 3) = 465 MeV, and Mb = 1.1 GeV. 


IV. DYNAMICS OF COUPLED CHANNELS AND AN ADMIXTURE OF NEIGHBORING STATES 


In the CC approach X(3872) is treated as the original 2 3 Pi ( QQ ) state, shifted down by the CC interaction with 
the DD* {DD*) channel. Moreover, in Ref. [2(| explicit parameters of this interaction were found, which put the pole 
of X(3872) in the multi-channel Green’s function exactly at the position of the open channel threshold. It is also 
important that within our approach the parameters are calculated and found to lie in a narrow range, i.e., they are 
not fitted. Below we use these fixed parameters to predict the admixture generated by the CC interaction, to the 
original 2 3 Pi state. 

Leaving the details of the CC formalism to the Appendix and the original papers na-m, we recapitulate briefly 
the essence of the mixing problem. We define the extra (open)-channel interaction by the function Vcc (q,q',E) 
in momentum space, which depends on the total energy E and contains all thresholds of the open channels. Then 
we are writing the full Green’s function G^q \ which describes the evolution of a QQ pair (produced in the e + e~ 
annihilation or the B decay process), interacting in the intermediate states with DD* {DD*) and finally ending up 
again as a QQ pair (in a similar e + e _ annihilation or the B decay process), as 


r (BB) 

U QQ 


n,m 


E-E 




( 22 ) 


where the matrices are defined as E = E n 6 nm , and u> = w nm with 

Wnm(E) = j i>i{q)V cc {q,q',E)ip m (q ')( 23 ) 

Here the symbol B in ( Bip n ) stands for an operator, defining the production (annihilation) process. In the case when 
the pure QQ state is already formed, one can put B = 1 , while in the e + e _ annihilation process the operator B is 
proportional to 7 M . 

Note, that G^q' 1 contains all CC contributions in the intermediate states, but the experimental conditions require 

these contributions to be projected onto the QQ states. We now include in G^q 1 the possibility to emit any particle 
( 7 , 7 r,...) in the intermediate state. As shown in Ref. |2l| (see Appendix), this can be done inserting the self-energy 
part with the photon loop to E and w: E + w —> E{ 7 ) T- w( 7 ), so that the probability of this emission is given by the 
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discontinuity of the self-energy part of Gqq 1 ' , e.g. 


Y- ( (E) = l A G, 


y(BB) 

2 ~"QQ 

= £(^») 

n,m,l 


1 


E — E + ui 


— A(E(7) + w(7)) m ; 


1 


E — E + w 




lq 


= ^ak^ky.A(E+ w)kni>na\i- 


k,n 


(24) 


In this way one can start in e + e , or the B decay, with the creation of any combination {P'0 n }, which finally, at the 
moment of the 7 detection, becomes the set of the states (ipk + 7 ) with the probability amplitude ak(Bip n ). 

The amplitude ak is given by 


Ofe 



1 

E — E + w 


{Btpi). 

kl 


(25) 


Since P is diagonal, P)*, = EiSik, all matrix structure is due to w. This matrix was calculated in Ref. plj for the 
1 ++ states, giving the values shown in the Table m in the Appendix, here one can see, that only two 1 ++ states are 
strongly connected by uu nm , namely l 3 Pi and 2 3 Pi, while higher states yield very small admixtures. 

Defining the amplitudes 01 and <22 for the 1 3 Pl and 2 3 Pi states, respectively, one has 

ar (E 2 + w 2 2-E)§fe-w 21 

„ ~ JV, | TT, Bl/)! ' 

02 Ei+wu - E- w 12 5 ^ 

With the use of the RSH we have E 2 = 3949 ± 10 MeV, Ei = 3523 ± 10 MeV and then from Table Hill 

E 2 + w 22 = (38 50 ± 10)MeV, Pi + tu n = (3203 ± 10)MeV. (27) 

From Eq. (l26l) one can see that the solution for ai/a 2 is sensitive to the original value of E 2 , but depends weakly 
on and we take for this ratio the same value as for the e + e _ process, namely |V’i( 0 )| 2 /|i/> 2 ( 0)| 2 ~ 0.85. 

Inserting these numbers into Eq. (l26l) and taking E 2 = 3949 MeV from Table U (the contribution from the D*D* 
channel is neglected) one obtains the value of cx at the resonance energy, Pr = 3872 MeV (P^i/ Bij) 2 = 0): 

cx(Pr) = —(P = 3872 MeV) = 0.183. (28) 

02 


However, this number can be considered as a lower limit, since we have not accounted for the contribution of the 
closed D*D* channel to the original levels Pi and E 2 . The estimate of these D*D* contributions to the 2 ++ ( 2 3 P 2 ) 
and 0 ++ (2 3 Po) states was made in Ref. [2(j for the same CC parameters as for X(3872). It yields a shift around —30 
MeV, in good agreement with the P(3930) as the 2 3 P 2 state. In this case in Eq. (l28l) E 2 = 3920 MeV and a larger 
value of cx is obtained, 


c x = — (P = 3872 MeV) « 0.215 
0-2 

As a result, we chose to accept the value of cx from the range, defined by Eqs. (051) and (051) 


c x = 0.200 ±0.015. 


(29) 


(30) 


V. THE COUPLED-CHANNELS CORRECTIONS TO THE X(3872) WAVE FUNCTION 

In the previous section the effects from open DD* channels were calculated, resulting in the contribution from the 
1 3 Pl state to the normalized w.f. of A(3872), which can be represented as 


¥>(*(3872)) = y^-4 ^(2 3 Pi) ±c x ^(1P 3 Pi) 


(31) 







As seen from Eq. (1301) the mixing parameter cx at the energy E r = 3.782 GeV lies in a narrow range; this admixture 
was not fitted, but calculated within the string breaking mechanism. 

Moreover, in the many-channels mechanism the sign of the continuum component cannot be arbitrary anymore and 
here we choose positive cx: it takes place if the radial w.f. P 2 p(r) is chosen to be negative at small distances (for 
r < 0.5 fm), while the radial w.f. of the ground state Rip(r) is positive in the whole region. Just for this choice a 
negative value for the overlap integral in Eqs. (HU El was obtained. 

Due to the CC admixture in the A(3872) w.f., the m.e. in the radiative decays of X (3872) can be presented as a 
superposition, 


(X(3872)\g(r)\J/il>) = y/l - c* Ar + c x /i, 

(A(3872)| ff (r)|0(25)) = \j\-c\ K 2 + cxh , (32) 

where I n and K n (n = 1,2) are defined in Eq. ([3 and the radial w.f. of "0(3686) is taken as R(2S)(r), i.e., without 
mixing with the ID state. 

For A'(3872) two different cc potentials are used below. In first case the w.f. of all states involved in the radiative 
decays of A(3872), are defined by the RSH with the linear+GE potential, Ve(r), given in Eq. ([9]). Then in the m.e. 
Eq. (1321) the two terms have different signs and therefore the magnitude of this m.e. is smaller as compared to the 
number given in Eq. (1151) . Due to this cancellation this m.e. appears to be very small, even for the minimal value of 
cx(min) = 0.185 from Eq. (1251) : 


|<X(3872,2 3 P 1 ))|, 7 (r)|J/0>| < 0.08 GeV" 1 , (33) 

giving a small partial width Ti = 3.0 keV. 

On the contrary, for the transition A(3872) —► 0(25)7 its m - e - increases to 3.30 GeV -1 and T 2 = 94 keV is 
large. Thus, within the CC approach, taking the linear+GE potential for A(3872), one obtains a small partial width 
Ti < 3 keV and very large value for the ratio 1.0. 

The situation is different, if we assume that the GE interaction is totally suppressed for higher charmonium states, 
i.e., the cc dynamics is determined by the confining potential. (Note that such a suppression does not exclude 
the use of the hyperfme interaction, having a very small range [28J, for the higher 5-wave charmonium states, like 
0(4040), concentrated near the origin). As a result, the m.e. of the radiative transitions, X(3872) —> J/0 7 and 
A(3872) —> 0(3686)7, are defined by the “mixed” m.e., in which the 2 P w.f. is calculated with the confining 
potential, while the w.f. of J/0 and 0(3686) are calculated with the linear+GE potential. 

These m.e. are denoted as, 


K\ (mix) = J drr 2 P 2 p(lin)|ff(r)|i?i S (lin + GE), 

AT 2 (mix) = J drr 2 P 2 p(lin)|g(r)|i? 2 s(lin + GE). (34) 

It appears that the magnitude of the m.e. Ad (mix) is larger compared to the number given in Eq. m and for the 
un-regularized w.f. (the solutions of the SSE in Eq. (fl2l) l the m.e. are 

50(mix) = —0.693 GeV -1 , AT 2 (mix) = 2.053 GeV -1 . (35) 

Instead, below we use the m.e. A,(mix) (i = 1,2), calculated for the regularized w.f. with the use of the Einbein 
approximation: 

Ad (mix) = -0.664 GeV -1 , AT 2 (mix) = 2.475 GeV -1 . (36) 

For the calculated Ad (mix) the partial width Ti strongly depends on the value of the mixing parameter c x , even if 
it is taken from the narrow range given in Eq. (l30ll . Using the m.e. I m (which refers to the IP state) from Tabic ITIll 
and AV(mix) from Eq. (1351) . one has 

Ti = T(AT(3872) J/0 7 ) = (37 ± 6 ) keV for c x = 0.200 ± 0.015. (37) 

Owing to the CC effects the value of Ti appears to be ~ 2.2 times smaller than in SCA. While changing the admissible 
value of cx from the minimal to the maximum values, i.e., by 15%, the partial width IJ increases from 31 keV to 43 
keV, i.e., 1.4 times. It means that a precise knowledge of Id from experiment would put a strong restriction on the 
mixing parameter cx at the resonance energy Ar = 3.872 GeV. 
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For the radiative decay X(3872) —> ^( 3686)7 we have used the radial w.f. i? 2 s(?') for t/’(3686). Then taking the 
m.e. I 2 S,ip from Table Hill and K 2 (mix) from Eq. (15C>1) . one obtains 

T 2 = (72 ± 1.0) keV for c x = 0.200 ± 0.015. (38) 

Its value weakly depends on the mixing parameter c x . However, this width may be smaller, if the 2 S — ID mixing 
is taken into account in the w.f. of ^>(3686) [H|. 

Thus in the CC approach with the minimal value of the mixing parameter, c x = 0.185, we have 

Ti = 42.7 keV, T 2 = 70.5 keV, R in = 1.65, (39) 

while for the maximal value, c x = 0.215, 

Ti = 30.8 keV, T 2 = 73.2 keV, R^ = 2.38, (40) 

i.e., this partial width Ti is by 40% smaller than for c x (min) = 0.185. The ratio of the widths also changes in a rather 
wide range: 


i?i/, 7 (th) = 2.02 ± 0.36, (41) 

being 2.4 times larger than that in the SCA Eq. (l20l) . This value is in a good agreement with the LHCb data given 
in Eq. ©. Also its lower limit (for c x = 0.185) does not contradict the Belle result, < 2.1. 

Thus, in our analysis the agreement with the LHCb result was obtained under two conditions: when the GE 
interaction is suppressed and the CC admixture to the w.f. of X(3872) is not small, ~ 20%. For the partial width Ti 
we predict the value 37 ± 6 keV. 


VI. CONCLUSIONS 

The cc dynamics in X(3872) with J PC = 1 ++ was studied within the CC approach, where a coupling to the DD* 
channels is determined by the string-breaking mechanism and the mixing parameter c x (E l = Er) = 0.200 ±0.015 was 
calculated. Two different variants of the cc interaction were considered for the higher charmonium states. 

1. The linear±GE potential used is the same for the lower and the higher charmonium states. In this case the 
partial width Ti = r(X(3872) —> J/i/j'y) ~ 3.0 keV is small and gives rise to a very large value of the ratio 
R^ 1.0 in the CC approach. On the contrary, in SCA the ratio R ^ 7 = 0.84 is not large. 

2. The GE interaction is supposed to be suppressed for higher states. For this dynamics the charmonium spectrum 
above the DD threshold can be described with a good accuracy, if perturbative corrections in the c-quark mass 
are neglected. 

If the suppressed GE potential is used, then the partial width Ti increases, being very sensitive to the value of the 
mixing parameter c x . For c x = 0.215 we find Ti = 30.8 keV, T 2 = 73.2 keV, and the calculated ratio of the partial 
widths R,/, 7 (max) = 2.38 is in good agreement with the LHCb result. The predicted minimal value, i?^ 7 (min) = 1.65, 
does not contradict the Belle measurements with R^, 7 (Belle) <2.1. 
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Appendix A: The coupled-channel mechanism 

Here we essentially use the formalism, suggested in Ref. [44] and developed further in Refs. [20j, [2lJ, and below we 
partly reiterate the material from Ref. [21] for the convenience of the reader. 

We use here the string-decay Lagrangian of the 3 Pq type for the decay cc -A ( cq){cq ) [20[ : 


£sd = 


d X , 


(Al) 
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where the light quark bispinors are treated in the limit of large m c , which allows us to go over to the reduced (2 x 2 ) 
form of the decay matrix elements. Moreover, to simplify the calculations the actual w.f. of the cc states, calculated 
in Ref. [21] with the use of the RSH, are fitted here by five (or three) oscillator w.f. (HO), while the D meson w.f. 
is described by a single HO term, which provides a few percent accuracy with the parameter /?2 — 0.48. In this case 
the factor M u in JAD is calculated to be M u = A\/2a / -Jtt /3 2 — 0.8 GeV (see below), which produces the correct total 
width of ^(3770). The transition m.e. for the decays (cc) n —> {D n2 D n3 ) are denoted here as n —> ri 2 ,n 3 , and in the 
2 x 2 formalism this m.e. reduces to 


J n 


, 0 ) = 


VnI. 


d 3 q 

( P + q ^ n2 


(A2) 


The factor 2/123 contains a trace of the spin-angular variables (for details see Refs. [2(1 [Il[), and 7 = + u^v D +e 0 ) ’ 

where the average over the Dirac denominator contains the light quark mass m q , the scalar Us = err and the vector 
Vb(r) = — potentials, and Cqi which is the Dirac eigenvalue in the static potential, created by the heavy quark: 


a-p + 13 (m q + C/ s (r))V>(r) = (e 0 - Vb(r))^(r) . 


(A3) 


Knowing the string tension a = 0.18 GeV 2 and the averaged momentum distribution 02 in the w.f. of a heavy-light 
meson, the value M u = was calculated in Ref. (2Qj. For the 15, 2 S, and 3 S charmonium states M u is different 
and equal to 0.65, 0.80, 1.10, respectively. Using the averaged values of the scalar and vector potentials, one obtains 
7 Ri 1.4. 

The intermediate decay channels, like DD*, induce an additional interaction “potential” Vcc(Q: Q 1 , E) (the quota¬ 
tion marks imply nonlocality and energy dependence of this potential): 


Vcc(q,q',E) 
Xn 2 n 3 (<b P) 



d 3 p X n2n3 (q,p)Xl an3 (q',p) 
(2n) 3 E-E n2n3 (p) 
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yi23(q,p)lpn 2 (q - p)^n 3 {q - p)- 


The new Hamiltonian with account of the CC interaction is 


(A4) 

(A5) 


H = H 0 + Vcc (A 6 ) 

where Hq is given in Eqs. m - At this point it is important to stress that Vcc contains all thresholds of open 
channels, and to treat the spectrum of the Hamiltonian H rigorously, one should exploit the Weinberg eigenvalue 
formalism, developed for this purpose in Ref. [2Cj, second reference. However, for the states below all thresholds (or 
neglecting the open channel amplitudes in first approximation) one can use an expansion in the complete set of the 
eigenvalues of Hq. One can show, as in Ref. [20[; that the exact Weinberg formalism reduces to this expansion for the 
real eigenvalues at or below thresholds. Using that, one obtains for the Green function an expansion, 


Gqq{ 1,2 ;E) = J2 MV(E -E + w)~^ m (2), 

n,m 


(A7) 


where the matrix element of w is denoted by w nm and given by 


^nm {E) 


d 3 q d 3 q' 
(27t) 3 (27t) 3 


^n(q)V C c{q,q' ,E)^ m (q') 


d 3 p \ ' dnn 2 n 3 {p)Jmn 2 nAp) 


{2nf 


E 


E — En 


,(P) 


and the energy eigenvalues are obtained from the determinant condition 

det(£’ — E — w) — 0, 


(AS) 


(A9) 


where (E") nm — E n 8 nrn 
Defining as in Eq. 01 the effective QQ w.f., 




k,l 


1 


E — E + w 


(Bipi) = y^afc'i/’fc, 

kl u 


(A10) 
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one has to calculate the resulting amplitudes a-k in the expansion of X(3872) in terms of the QQ states. The values 
of w nm are given in Table ITVl and one can see that only two states, |1P) and |2P), are important. 

Keeping only two QQ eigenfunctions, one has 


ai 


a 2 


E 2 


- E + w 2 2 
det 


(P^i) 


W21 

det 


(Bih) 


Ei 


- E + w\i 
det 


{Bih) - 1 ), 


where det = det(E — E + w). Then one obtains the ratio 


01 

0-2 


—W 21 + (E 2 + W 22 — E) 


W 1 ) 

(B^a) 


Ei + w 11 — E — w 12 


(Bbi) 

(Bb 2) 


(All) 

(A12) 


(A13) 


We calculate this ratio for the energy Pr = 3872 MeV, taking the values of Wik from the Table HV1 Note, that we have 
corrected the sign of W 12 in accordance with the relative signs of the wave functions |1P) and |2P). The resulting 
values of cx = ^-(3872 MeV) are given in Eqs. (I28M30I) . 

From Table m one can see that the 3 3 Pi state gives a negligible admixture. The values of w mn for E = Er are 
computed according to Eq. (IA8I) with the w.f., calculated in Ref. [20], and are given in Table ITVl 


TABLE IV: The m.e. w nm (in GeV) between n 3 Pi and m 3 Pi states for two approximations of the exact w.f. 


n m 

11 

1 2 

2 2 

3 2 

Wnm for 5 HO 

-0.320 

0.122 

- 0.099 

-0.0003 

Wnm for 3 HO 

-0.319 

0.121 

- 0.098 

-0.0011 
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